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The Weakly Interacting Massive Particle (WIMP) is the main candidate for the relic dark matter. 
A set of exclusion curves currently obtained for cross sections of the spin-dependent WIMP-proton 
and WIMP-neutron interaction is given. A two-orders-of-magnitude improvement of the sensitivity 
of the dark matter experiment is needed to reach the SUSY predictions for relic neutralinos. It 
is noted that near- future experiments with the high-spin isotope 73 Ge can yield a new important 
constraint on the neutralino-neutron effective coupling and the SUSY parameter space. 

A. Introduction 

Nowadays the main efforts in the direct dark matter search experiments are concentrated in the field of 
so-called spin-independent (or scalar) interaction of a dark matter particle or the Weakly Interacting Massive 
Particle (WIMP), with a nucleus. The lightest supersymmetric (SUSY) partilce (LSP) neutralino is assumed 
here as a best WIMP candidate. It is believed that this spin-independent (SI) interaction of dark matter 
(DM) particles with nuclei makes a dominant contribution to the expected event rate of detection of these 
particles. The reason is the strong (proportional to the squared mass of the target nucleus) enhancement of 
SI WIMP-nucleus interaction. The results currently obtained in the field are usually presented in the form of 
exclusion curves (see for example Fig.^). For the fixed mass of the WIMP the values of the cross section due 
to scalar elastic WIMP-nucleon interaction located above these curves are already excluded experimentally. 
There is also the DAMA closed contour which corresponds to the first claim for evidence for the dark matter 
signal due to the positive annual modulation effect 

In the paper we consider some aspects of the spin-dependent (or axial-vector) interaction of the DM 
WIMP with nuclei. There are at least three reasons to think that this spin-dependent (SD) interaction could 
also be very important. First, contrary to the only one constraint for SUSY models available from scalar 
WIMP-nucleus interaction, the spin WIMP-nucleus interaction supplies us with two such constraints (see for 
example and formulas below). Second, one can notice 0,0 that even with a very accurate DM detector 
(say, with sensitivity 10 -5 events/day/kg) which is sensitive only to the WIMP-nucleus scalar interaction 
(with spinless target nuclei) one can, in principle, miss a DM signal. To safely avoid such a situation one 
should have a spin-sensitive DM detector, i.e. a detector with spin-non-zero target nuclei. Finally, there is 
a complicated (and theoretically very interesting) nucleus spin structure, which possesses the so-called long 
q-tail form-factor behavior for heavy targets and heavy WIMP. Therefore, the SD efficiency to detect a DM 
signal is much higher than the SI efficiency, especially for the heavy target nucleus and WIMP masses 

B. Zero Momentum Transfer 

A dark matter event is elastic scattering of a relic neutralino \ ( or x) from a target nucleus A producing 
a nuclear recoil Er which can be detected by a suitable detector. The differential event rate in respect 
to the recoil energy is the subject of experimental measurements. The rate depends on the distribution of 
the relic neutralinos in the solar vicinity f(v) and the cross section of neutralino- nucleus elastic scattering 
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FIG. 1: WIMP-nucleon cross section limits in pb for scalar (spin-independent) interactions as a function of the WIMP 
mass in GeV. Shown are contour lines of the present experimental limits (solid lines) and of projected experiments 
(dashed lines). Also shown is the region of evidence published by DAMA. The theoretical expectations are shown by 
scatter plots (circles and triangles are from 0,0) an d by grey closed region 3]. From Q. 



14L 11.1 lid . The differential event rate per unit mass of the target material has the form 
dR 



dE R 



N £x f maX dvf(v)v^(v,q 2 ). 



(1) 



The nuclear recoil energy = q 2 /(2Ma) is typically about 10 _6 m x and N = JV/A is the number density 
of target nuclei, where J\f is the Avogadro number and A is the atomic mass of the nuclei with mass Ma- 
The neutralino-nucleus elastic scattering cross section for spin-non-zero (J ^ 0) nuclei contains coherent 
(spin- independent, or SI) and axial (spin-dependent, or SD) terms [^ IrRllsj : 



da J 
~dq~ 2 
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The normalized non-zero-momentum-transfer nuclear form-factors 

F^sM 2 ) = (*&,si(0) = 1), (3) 



? 2 i„2\ _ ^SD,SI V^ J 1^2 

are defined via nuclear structure functions |8L 1171 ll 



siM = E l(^IICL(g)||J)| 2 ^|(J||Co( g )||J>| 2 , (4) 

L even 

^SD(g) = E {\{N\\T£ m {q)\\N)\ 2 + \{N\\Cl{q)\\N)\ 2 ). (5) 



Lodd 



The transverse electric T el5 (q) and longitudinal C 5 (q) multipole projections of the axial vector current 
operator, scalar function Cl{m) are given in the form 



Cdq) = E c o3L(qri)Y L (fi), C (q) = ^2c j (qr i )Y (f i ), 



clcons 



where a 01 = a n ± a p (see (JTHJ) and M L . LI {qri) = 3L'{qri)[Y L '{f l )ai} L lalHHlal- The nuclear SD and SI 
cross sections at q = (in J2J) have the forms 



^ ,0) - (2J+1) =^ Af7sl(0) ' ,(n 
,2 c /T^ /i ,,2 



V^(0)<S£> +sign(a p a n )y / a s " D (0)(S^) I . (8) 



A fr\\ MS SD (0) 4^(J+l)r A A 

= f A (J + l) 

Here /i^ = — — is the reduced neutralino- nucleus mass. The zero- momentum-transfer proton and 

m x + Ma 

neutron SI and SD cross sections 

4M = ^4, Co ^ c r^E c «4 p ' n) ; ( 9 ) 

aim = l2 ^ a 2 n fl „ = ^ M W a P =E^W (10) 

depend on the effective neutralino-quark scalar C q and axial- vector _4 g couplings from the effective Lagrangian 

£cff = E (A ' ^7m75X ' 97 M 759 + C g • XX ' 5?) + ••• (H) 

and on the spin (A^< n) ) and mass {fq P ' n ^) structure of nucleons. The factors Aq P ' n ^ parametrize the quark 
spin content of the nucleon and are defined by the relation 2A^ ra ' p ^s p = (p, s|^ g 7 M 75^ 9 \p, s)( Pi „). The total 
nuclear spin (proton, neutron) operator is defined as follows 



where i runs over all nucleons. Further the convention is used that all angular momentum operators are 
evaluated in their z-projection in the maximal Mj state, e.g. 

(S) = (N\S\N) = (J,Mj = J\S z \J,Mj = J). (13) 

Therefore (S p / n )) is the spin of the proton (neutron) averaged over all nucleons in the nucleus A. The cross 
sections at zero momentum transfer show strong dependence on the nuclear structure of the ground state 

000. 

The relic neutralinos in the halo of our Galaxy have a mean velocity of (v) ~ 300 km/s = 10~ 3 c. When the 
product (7 m ax-R <C 1, where R is the nuclear radius and <7 max = Sf^AV is the maximum momentum transfer in 
the x^4 scattering, the matrix element for the spin-dependent \A scattering reduces to a very simple form 
{zero momentum transfer limit) [2(1 l2l|: 

M = C(N\a p S p + a n S n \N) ■ = CA(N\3\N) ■ s*. (14) 

Here s x is the spin of the neutralino, and 

(N\a p S p + a n S n \N) (N\(a p S p + a n S n ) ■ JjjV) 

(N\J\N) J(J + 1) ' [ j 

It is seen that the \ couples to the spin carried by the protons and the neutrons. The normalization C 
involves the coupling constants, masses of the exchanged bosons and various LSP mixing parameters that 
have no effect upon the nuclear matrix element 22]. In the limit of zero momentum transfer q = the spin 
structure function 10) reduces to 

Ssd(0) = — — A 2 J(J + 1). (16) 

TT 

Perhaps the first model to estimate the spin content in the nucleus for the dark matter search was the 
independent single-particle shell model (ISPSM) used originally by Goodman and Witten 23] and later 
in [llj, I24L |25| . The ground state value of the nuclear total spin J can be described by those of one extra 
nuclcon interacting with the effective potential of the nuclear core. 

There are nuclear structure calculations (including non-zero-momentum approximation) for spin- 
dependent neutralino interaction with helium 3 He 0; fluorine 19 F 000;sodium 23 Na 00 26 , 27 1; 
aluminium 27 A 10 ; silicon 29 Si 0,00; chlorine 35 C1 0; potassium 39 K 0; germanium 73 Ge I3Es|; 
niobium 93 Nd Ql iodine 127 I Q; xenon 129 Xe and 131 Xe 0,00; tellurium 123 Te 0and 125 Te 
0]; lead 2O8 Pb0,0]. The zero-momentum case is also investigated for Cd, Cs, Ba and La in0,0,0]. 

There are several approaches (advocated by a few groups of researchers) to the more accurate calculation 
of the nuclear structure effects relevant to the dark matter detection. To the best of our knowledge an almost 
full list of the models includes the Odd Group Model (OGM) of Engel and Vogel 0] and their extended 
OGM (EOGM) 00; Interacting Boson Fermion Model (IBFM) of Iachelk^Krauss, and Maino 0; 
Theory of Finite Fermi Systems (TFFS) of Nikolaev and Klapdor-Kleingrothaus j3jj; Quasi Tamm-Dancoff 
App roximation (QTDA) of Engel different shell model treatments (SM) by Pacheco and Strottman 



, by Engel, Pittel, Ormand and Vogel |2jj and Engel, Ressell, Towner and Ormand, [2j, by Ressell et 
al. and Ressell and Dean 0; by Kosmas, Vergados et al. {lH Efil Isij : so-called "hybrid" model of 
Dimitrov, Engel and Pittel 0] and perturbation theory (PT) based on calculations of Engel et al. . 



C. Spin constraints 



For the spin-zero nuclear target the experimentally measured event rate Q of direct DM particle detection 
via formula is connected with zero-momentum WIMP-proton (for the neutron the cross section is the 



same) cross section (jTJ). The zero momentum scalar WIMP-proton (neutron) cross section <7gj(0) can be 
expressed through effective neutralino-quark couplings C q Ijf f(l by means of expression J^JJ • These couplings 
C q (as well as A q ) can be directly connected with the fundamental parameters of a SUSY model such as 
tan/3, M^a, Mi masses of sfermions and Higgs bosons, etc. Therefore experimental limitations on the spin- 
independent neutralino-nucleon cross section supply us with a constraint on the fundamental parameters of 
an underlying SUSY model. 
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FIG. 2: Full set of currently available exclusion curves for spin-de pen dent WIMP-proton cross sections (a ^ as a 



function of WIMP mass) . The curves are obtained from |3£ 
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In the case of the spin-dependent WIMP-nucleus interaction from measured differential rate QJ one first 
extracts limitation for cr<^ D (0) and therefore has in principle two constraints [a] for the neutralino-proton 
a p and neutralino- neutron a n spin effective couplings as follows from relation Q. From Q) one can also 
see that contrary to spin-independent case JSJl there is no factorization of the nuclear structure for cr^ D (0). 
Both proton (S^ ) and neutron (S^) spin contributions simultaneously entering formula for the SD 
WIMP-nuclear cross section cr^ D (0)- 

In the earlier considerations based on the OGM Q, Q one assumed that the nuclear spin is carried by 



WIMP-neutron axial coupling limits 

x> 




FIG. 3: Full set of currently available exclusion curves for spin-dependent WIMP-neutron cross sections (o"sd versus 
WIMP mass). The curves are obtained from 0, ED, El EE El Q Note that the NAIAD curve corresponds to the 
sub-dominant contribution, extracted from the p-odd nucleus Na (compare with the relevant NAIAD curve in Fig.[5J 



the "odd" unpaired g roup of protons or neutrons and only one of either (S„) or (S^ ) is non-zero (the same 
is true in the ISPSM [111 l23l 124 , l25|). In this case all possible target nuclei can naturally be classified into 
n-odd and p-odd groups. The current experimental situation for the spin-dependent WIMP-proton cross 
section is given in Fig. [21 The data are taken from experiments BRS, (Nal, 1992) 0E1, BPRS JCaFo, 

1996) 0, DAMA (NAI, 1996) 0, DAMA (CaF 2 , 1999)0, 0, 
ELEGANT (CaF 2 , 1998) 0, ELEGANT (Nal, 1999) HQ, Tokio 



1993) |3J, EDELWEISS (sapphire 



UKDMS (Nal, 1996) UAW, 
(LiF, 1999, 2002) 0, El E 

(C„F m , 2002) 0, ANAIS (Nal, 2002) |5^| and NAIAD (Nal, 2003) |5||. The current experimental situation 



53j,|5|, SIMPLE (C 2 C1F 5 , 2001) |55|, CRESST (A1 2 3 , 2002) |56j, PICASSO 



for the spin-dependent WIMP-neutron cross section is given in Fig. |3| The data are taken from the first 



experiments with natural Ge (1988, 1991) |60]|6l(, xenon (DAMA/Xe-0,2) |62j, |63j, |6J and sodium iodide 
(NAIAD) 0. In the future one can also expect exclusion curves for the SD cross section, for example, from 
the EDELWEISS and CDMS experiments with natural germanium bolometric detectors. 

From Fig. 0] one can conclude that an about two-orders-of-magnitude improvement of the current DM 
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FIG. 4: The same as in Fig. [5] but with the theoretical scatter plot from obtained in the effMSSM with all 
coannihilation channels included (green circles) and with 0.1 < Qh 2 < 0.3 (black triangles). The triangle-like 
shaded area is taken from 31. 



experiment sensitivities is needed to reach the SUSY predictions for the crg D provided the SUSY lightest 
neutralino (LSP) is the best WIMP particle candidate. There is the same situation for the <7grj- 

Further more accurate calculations ofj!E3,E10E10EiEiElEiEi demonstrate that contrary 
to the simplified odd-group approach both (S^ ) and (S^) differ from zero, but nevertheless one of these 
spin quantities always dominates ((S^) <C (S^), or (S^) <C (S^)). If together with the dominance like 
(^p(n)) ^ (^«(p)) one wou ld have the WIMP-proton and WIMP-neutron couplings of the same order of 
magnitude (not a n M <C GW n )), the situation could look like that in the odd-group model. Nevertheless it 
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FIG. 5: The DAMA/Nal region from the WIMP annual modulation signature (57986 kg day) in the (£asi, £ct S d) 
space for 40 GeV< m WIMP <110 GeV covers all four particular couplings (9 — 0, 9 — 7r/4, 9 — tt/2 and 9 — 2.435 
rad) reported in Scatter plots give correlations between o-g : (0) and <r SD in effMSSM (£ = 1 is assumed) with 

the requirement for neutralino relic density 0.1 < Q.h 2 < 0.3 and all coannihilation channels included J2J. Red stars 
correspond to an assumption that m LSP < 150 GeV. 

was shown in J35] that in the general SUSY model one can meet a case when a n ( p ) <C a p ( n )- To solve the 
problem (to se par ate SD proton and neutron constraints) at least two new approaches were proposed. As 
the authors of 65] claimed, their method has the advantage that the limits on individual WIMP-proton and 
WIMP-neutron SD cross sections for a given WIMP mass can be combined to give a model-independent 
limit on the properties of WIMP scattering from both protons and neutrons in the target nucleus. The 
method relies on the assumption that the WIMP-nuclear SD cross section can be presented in the form 

V 

t J , where tfg D U and o"sdU are auxiliary quantities, not directly connected 



^sd(O) 



r SDlyl 



± 



SDlA 



with measurements. Furthermore, to extract, for example, a constraint on the sub-dominant WIMP-proton 
spin contribution one should assume the proton contribution dominance for a nucleus whose spin is almost 
completely determined by the neutron-odd group. It may seem useless, especially because these sub-dominant 
constraints are always much weaker than the relevant constraints obtained directly with a proton-odd group 
target. 

Another approach of |66l is based on introduction of another auxiliary quantity asr> = 12— (a? + a?), 
where tan 9 = a n /a p . With these definitions the SD WIMP-proton and WIMP-neutron cross section can be 
given in the form Og^(0) = ctsd sin 2 #(cos 2 9). In Fig.[S]the WIMP-nucleon spin-mixed and scalar couplings 
allowed by annual modulation signature from the 100-kg DAMA/Nal experiment (57986 kg day) are given 
by filled region. The region (at 3 a C.L.) in the (£ctsi, £ctsd) s P ace f° r 40 GeV< m wlMP <110 GeV covers 
all four particular couplings [0 = 0, 6 = n/4, 9 = n/2 and 9 = 2.435 rad) reported in 6.6]. Scatter plots 
give fTgj(O) versus osd in effMSSM with 0.1 < flh 2 < 0.3 and all coannihilation channels included from Q 
under the assumption of £ = 1. Red stars correspond to the same calculations but with m LSP < 150 GeV. 
In this mixed case the limits for the spin couplings depend on assumptions about the scalar coupling, and 
the relevant exclusion curve for the spin-dependent WIMP-proton cross section (not given in Fig. can n ot 
be simply extracted from these mixed results of j6?j |. 



D. The role of the germanium-73 isotope 



Comparing the number of exclusion curves in Fig. |2 and Fig. [31 one can easily see that there are many 
measurements with p-odd nuclei and there is a lack of data for n-odd nuclei, i.e. for crg D . Therefore 
measurements with n-odd nuclei are needed. From our point of view this lack of cjg D measurements can 
be filled with new data expected from the HDMS experiment with the high-spin isotope 73 Ge [f^. This 
isotope looks with a good accuracy like an almost pure n-odd group nucleus with (S n ) S> (S p ) (Table [I]). 
The variation of the (S p ) and (S„) in the table reflects the level of inaccuracy and complexity of the current 
nuclear structure calculations. 



TABLE I: Zero-momentum spin structure (and predicted magnetic moments n) of the 73 Ge nucleus in different nuclear 
models. The experimental value of the magnetic moment given in the brackets is used as input in the calculations. 



73 Ge (L., = G 9/2 ) 


(Sp> 


(Sn) 


fi (in /ijv) 


ISPSM, Ellis-Flores [24,££] 





0.5 


-1.913 


OGM, Engel-Vogel [34J 





0.23 


(-0.879)cx P 


IBFM, Iachello at al. [22] and UU 


-0.009 


0.469 


-1.785 


IBFM (quenched), Iachello at al. [33] and [17] 


-0.005 


0.245 


(-0.879) cxp 


TFFS, Nikolaev-Klapdor-Kleingrothaus, [2§] 





0.34 




SM (small), Ressell at al. [1JJ 


0.005 


0.496 


-1.468 


SM (large), Ressell at al. [1JJ 


0.011 


0.468 


-1.239 


SM (large, quenched), Ressell at al. ^JJ 


0.009 


0.372 


(-0.879) OX p 


"Hybrid" SM, Dimitrov at al. [28] 


0.030 


0.378 


-0.920 



On the other hand, Fig. shows that for the ratio of a n to a p one can have the bounds 



0.55 < 



On 

n,, 



< 0. 



The scatter plots in Fig. were obtained in efTMSSM [2( when all coannihilation channels were included. 
The blue squares (black points) were calcualted with (without) the relic neutralino density constraint 0.1 < 
flh 2 < 0.3. Therefore in the model the couplings are almost the same and one can safely neglect the (S^)-spin 
contribution in the analysis of the DM data with the 73 Ge target (for which (S^ ) <C (S^)). 
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FIG. 6: The scatter plot (dots) gives the ratio of the neutralino-neutron a n and neutralino-proton a p spin couplings 
m the effMSSM 0. Boxes correspond to the relic neutralino density constraint 0.1 < Q.h 2 < 0.3 in the same model. 



We would like to advocate the old odd-group-like approach for experiments with germanium detectors. 
Of course, from measurements with 73 Ge one can extract not only the dominant constraint for WIMP- 
nucleon coupling a n (or crg r>) b ut also the constraint for the sub-dominant WIMP-proton coupling a p (or 
°sd) usm g t ne approach of |65| . Nevertheless, the latter constraint will be much weaker in comparison with 
the constraints from p-odd group nuclear targets, like 19 F or Nal. The fact illustrates the NAIAD (Nal, 
2003) curve in Fig. [21 which corresponds to the sub-dominant WIMP-neutron spin contribution extracted 
from the p-odd nucleus Na. 



E. Finite Momentum Transfer 

As becomes larger, the product qR ceases to be negligible and the hnite momentum transfer limit must 
be considered for heavier nuclei. With the isoscalar coupling constant ao = a n + a p and the corresponding 
isovector coupling constant a\ = a n — a n one splits S^ D (q) into a pure isoscalar, Sqo, a pure isovector, Sn, 
and an interference term, Sqi [13, |2fJ : 

Sin(q) = a 2 Q S*(q) + a^q) + ao^S^q). (17) 



The differential SD event rate has the form 



dq 2 



vdvf(v) 



8G^ 



(2J+ l)v 



(18) 



Comparing the differential rate l|18[> together with the spin structure functions of (|17() with the observed 
recoil spectra for different targets (Ge, Xe, F, Nal, etc) one can directly and simultaneously restrict both 
isoscalar and isovector neutralino-nucleon effective couplings an.i- These constraints will impose most model- 
independent restrictions on the MSSM parameter space free from any assumption of Perhaps, it 
would be the best to fit all data directly 1651 in terms of neutralino proton and neutron effective spin couplings 



l 0,l 



or a n (see for example analysis of [54^ and not to use such spin quantities as Cg^ and a 



'sd- 



Another attractive feature of the spin-dependent WIMP-nucleus interaction is the q-dependence of SD 
structure function 1)17(1 . One knows that the ratio of SD to SI rate in the 73 Ge detector grows with the WIMP 
mass 0,E|- The growth is much greater for heavy target isotopes like xenon. The reason is the different 
behavior of the spin and scalar structure functions with increasing momentum transfer. For example, the 
xenon SI structure function vanishes for q 2 ss 0.02 GeV, but the SD structure function is a non-zero constant 
in the region (Fig. [7J) . 
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FIG. 7: The 131 Xe structure function for a pure bino. The single-particle structure function has been normalized to 
5(0) = 1. From 0. 

As noted by Engel in the relatively long tail of the SD structure function is caused by nucleons near 
the Fermi surface, which do the bulk of the scattering. The core nucleons, which dominate the SI nuclear 
coupling, contribute much less at large q. Therefore the SD efficiency for detection of a DM signal is higher 
than the SI efficiency, especially for very heavy neutralinos. 



F. Conclusion 



The idea of this review paper is to attract attention to the role of the spin-dependent WIMP-nucleus 
interaction in the dark matter search experiments. The importance of this interaction is discussed. The 
fullest possible set of currently available exclusion curves for spin-dependent WIMP-proton and WIMP- 
neutron cross sections is given in Fig. Eland Fig. EI Nowadays about two-orders-of-magnitude improvement 
of the current DM experiment sensitivities is needed to reach the SUSY predictions for the Cg^ . It is noted 
that a near- future experiment like HDMS j^] with the high-spin isotope 73 Ge being an almost pure n-odd 
nucleus can fill in this gap and can supply us with new important constraints for SUSY models. 
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